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Abstract. The luminous infrared-loud quasar 
IRAS 13349+2438 was observed with the XMM- 
Newton Observatory as part of the Performance 
Verification program. The spectrum obtained by the 
Reflection Grating Spectrometer (RGS) exhibits broad 
(FWHM ~ 1400 km s” 1 ) absorption lines from highly 
ionized elements including hydrogen- and helium-like 
carbon, nitrogen, oxygen, and neon, and several iron 
L-shell ions (Fexvil - xx). Also shown in the spectrum 
is the first astrophysical detection of a broad absorption 
feature around A=16-17A identified as an unresolved 
transition array (UTA) of 2p - 3d inner-shell absorption 
by iron M-shell ions in a much cooler medium; a feature 
that might be misidentified as an Ovn edge when 
observed with moderate resolution spectrometers. No 
absorption edges are clearly detected in the spectrum. We 
demonstrate that the RGS spectrum of IRAS 13349+2438 
exhibits absorption lines from two distinct regions, one 
of which is tentatively associated with the medium that 
produces the optical/UV reddening. 
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1. Introduction 

IRAS 13349+2438 is an archetypal highly-polarized radio- 
quiet quasar at a redshift of z — 0.10764 (Kim et al. 1995). 
Since its identification as an infrared-luminous quasar 
(Beichman et al. 1986), this source has been extensively 
studied in the optical, infrared, and X-ray bands. In a 

* Based on observations obtained with XMM-Newton , an 
ESA science mission with instruments and contributions di- 
rectly funded by ESA Member States and the USA (NASA). 


detailed investigation of the optical and infrared spec- 
tra and polarization measurements, Wills et al. (1992) 
demonstrated that the nuclear spectrum exhibits two dis- 
tinct components; a highly-reddened component and a 
highly-polarized component that suffers much lower ex- 
tinction. Based on these observational facts, Wills (1992) 
constructed a simple and elegant picture of the nuclear re- 
gion of IRAS 13349+2438 in which the direct AGN radia- 
tion is attenutated through a thick dusty torus, while the 
observed highly-polarized light is produced by scattering 
in an extended bipolar region, either by warm electrons or 
by small dust grains. 

IRAS 13349+2438 was detected in the ROSAT All- 
Sky-Survey (Walter & Fink 1993; Brinkmann & Siebert 
1994), and has been the target of extended pointed obser- 
vations with ROSAT (Brandt, Fabian, & Pounds 1996), 
and with ASCA (Brinkmann et al. 1996; Brandt et al. 
1997). In general, the analyses of these data sets have 
suggested the presence of absorption edge features due 
to ionized species of oxygen (0 VIII and O vn), although, 
because of the low spectral-resolving-power capabilities 
of the available detectors on these missions, the precise 
nature of the soft X-ray spectrum has remained contro- 
versial. Siebert & Brinkmann (1999) self-consistently ac- 
counted for the effects of dust embedded in the warm ab- 
sorbing medium, and concluded that single zone models, 
both with and without internal dust, do not provide ad- 
equate fits to the combined, ROSAT , ASCA, and optical 
data sets. In particular, they find that a dust-free warm 
absorber model formally gives the best fit to the X-ray 
data, and conclude that the X-ray absorption and optical 
reddening must arise in spatially distinct regions. 

In this Letter, we present results from the first high- 
resolution X-ray observation of IRAS 13349+2438 with 
the XMM-Newton Observatory. The spectrum obtained 
with the Reflection Grating Spectrometer (RGS) shows 
a wealth of discrete spectral features, including the first 
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astrophysical detection of inner-shell 2 p - 3d absorption 
by M-shell iron ions in the form of an unresolved tran- 
sition array (UTA). From a detailed analysis of the rich 
absorption spectrum, we measure the column density and 
velocity field of the line-of-sight material. We show that 
the spectrum contains absorption features from regions 
with two distinct levels of ionization. The column density 
of the lower ionization component is consistent with the 
observed optical reddening, and we tentatively associate 
this component with the dusty torus. 

2. Observation and Data Reduction 

IRAS 13349+2438 was observed with the XMM-Newton 
observatory (Jansen et al. 2000) on 19 - 20 June, 2000 
during the Performance Verification phase for a total ex- 
posure time of 42 ks. The data obtained with the RGS 
(den Herder et al. 2000) were filtered through standard 
event-selection criteria using the XMM-Newton Science 
Analysis Software (SAS). The source spectrum was ex- 
tracted using both dispersion/cross-dispersion coordinates 
and dispersion/pulse-height coordinates, which include 
90% of the source photons. The background spectrum was 
generated using all of the events that lie outside the spatial 
mask. Wavelengths were then assigned to the dispersion 
coordinates. The current wavelength scale is accurate to 
within ~ 8 mA across the entire RGS band of A = 5 — 35 A 
(E = 0.35 -2.5 keV). 

The European Photon Imaging Camera (EPIC; Turner 
et al. 2000) MOS1 data were also processed with the SAS. 
The MOS2 detector was operated in FAST UNCOMP 
mode, for which reduction by the SAS is not possible as 
yet. Source events were extracted using a circular region 
of radius 45". A nearby source- free region with a radius 
of 3' was used for the background. 

Three modes of the Optical Monitor (OM; Mason et 
al. 2000) were used during the observation. The V grism 
(or optical grism; exposure time 3000 s), UV grism (expo- 
sure time 1000 s), and the UVW2 filter (~ 1500 — 3000 A; 
effective exposure time 6340 s). There is no significant 
variability in the UVW2 observations. Absolute flux cali- 
brations for the grism data were not available at the time 
of writing. The optical spectrum shows dear evidence of 
H/3 emission but the signal- to- noise is insufficient to de- 
termine the presence of any other emission lines (within 
an observer frame wavelength range of ~ 3000 — 6000 A). 
The UV spectrum (coverage ~ 2000 — 3500 A) is too weak 
for spectral extraction. 

3. Results of Spectral Fitting 

3.1 . Underlying Continuum Radiation 

To first obtain a rough characterization of the shape of 
the continuum, we use the EPIC-PN spectrum, which has 
the highest statistical quality and covers a broad range in 
energy. The 0.2 - 10 keV spectral region can be well-fitted 
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Fig. 1. The RGS first order spectrum of IRAS 13349+2438 
corrected for cosmological redshift (z = 0.10764). The error 
bars represent 1 <j Poisson fluctuations. 

with a phenomenological model that consists of a power- 
law and two black body components absorbed through a 
Galactic column density of iV^ aI = l.l x 10 2 ° cm -2 (Mur- 
phy et al. 1996). The EPIC-PN data require two black 
body components with temperatures of kT ~ 70 eV and 
~ 12 eV. The best-fit powerlaw photon index of T ~ 2.2 
is consistent with the value implied by the A SC A data 
(Brinkmann et al. 1996; Brandt et al. 1997). 

The MOS1 data from 0.3 to 10 keV, excluding the 
0.6 - 1.2 keV region where the RGS shows complex ab- 
sorption features, were best-fit using a single blackbody 
plus power-law (a \r 0 f 1*59 for 199 degrees of freedom), 
which is generally consistent with the PN data. The black- 
body temperature is kT ~ 1 10 eV. The best-fit photon in- 
dex is T = 2.1 with a column density of 1.4 x 10 20 cm -2 , 
which is only slightly higher than the Galactic value of 
1 x 10 20 cm" 2 . 

3.2. The Absorption Spectrum of IRAS 13349+2438 

The RGS spectrum shown in Figure 1 exhibits numer- 
ous absorption lines from a wide range in levels of ioniza- 
tion. The most prominent features in the spectrum are K- 
shell absorption lines of H- and He-like carbon, nitrogen, 
oxygen, and neon, and L-shell lines of Fexvil — XX. The 
spectrum also shows a broad absorption feature between 
A ~ 16 - 17 A (E ~ 730 - 770 eV). The observed location 
of the red “edge” of this feature is at A = 17.10 ± 0.05 A 
in the rest-frame of the quasar, and is close but undoubt- 
edly inconsistent with that of the O vil photoelectric edge 
(A = 16.78 A). The shape of the absorption trough to- 
wards the shorter wavelengths is also incompatible with 
that of a photoelectric edge. We identify this feature as a 
UTA of inner-shell 2 p - 3d resonance absorption lines in 
relatively cool, M-shell iron. The shape of this feature is 
strikingly similar to laboratory absorption measurements 
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of a heated iron foil (Chenais-Popovics et al. 2000), and 
also agrees well with our own calculations as described 
below. 

We adopt a continuum model similar to that inferred 
from the EPIC-PN data; i.e., the sum of a powerlaw and 
a blackbody component, the latter merely in order to 
parametrize empirically the soft X-ray spectral shape. We 
also fix the powerlaw photon index to T = 2.2 and the 
column density of cold material to the Galactic value. 
With the continuum model defined, we then apply ab- 
sorption components of H- and He-like C, N, O, and Ne, 
and Fe XVII - xxiv. Each ion is treated as a separate com- 
ponent in the spectral fit and contains all relevant reso- 
nance transtions from the ground-state and photoelectric 
edges in the RGS band. The line profiles are calculated 
accounting for thermal and turbulent velocity broadening. 
Transition wavelengths and oscillator strengths of the up- 
per levels were calculated with the atomic physics package 
HULLAC (Bar-Shalom et al. 1998), except for the wave- 
lengths of the strong Fe L resonance lines where we use 
laboratory measurement values described in Brown et al. 
(2000). We use photoionization cross sections from Verner 
et al. (1996). For the iron UTA absorption, photon impact 
processes involving L-shell to M-shell excitations in FeV 
- XVI are computed. All of the transitions 21 8 31 x - 21 7 
3jz+i ( x - i through 12) are taken into account, of which 
the 2p - 3d excitations are the most important. For the 
atomic structure, the most significant configuration mix- 
ings, which conserve the total angular momentum within 
the n = 3 shell (namely 3p 2 + 3s3d), are included. This 
approximation is expected to be adequate for analysing 
the presently observed unresolved absorption feature. A 
more detailed discussion of the UTA is presented in Be- 
har, Sako, & Kahn (2000). 

Each component is then fit for the ion column den- 
sity simultaneously with the black body continuum pa- 
rameters and the normalization of the powerlaw com- 
ponent. The best-fit black body temperature is kT ~ 
100 eV with a flux in the 5 - 35 A range of 4.0 x 
10“ 3 photons cm" 2 s~*. The flux in the powerlaw com- 
ponent is 2.4 x 10 -3 photons cm -2 s -1 in the same wave- 
length range. 

The observed widths of the absorption lines, as shown 
below, are much larger than both the instrument line 
spread function and the thermal broadening of a gas with 
kT ~ 10 eV, which is the expected temperature for a 
photoionized plasma at this level of ionization. With the 
current statistical quality of the RGS data, however, we 
are not able to constrain the turbulent velocities v tur b of 
the individual ion components. We, therefore, assume a 
uniform mean turbulent velocity field, keeping in mind 
that each ion can, in principle, exist in regions of differ- 
ent turbulent velocities. The derived ion column densities, 
therefore, may be uncertain to some degree, as quantified 
in the following section. 
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Fig. 2. Same as in Figure 1 with the best-fit model spectrum 
superimposed. 

We obtain a statistically acceptable fit to the RGS 
data with xl = 1*07 for 532 degrees of freedom. The con- 
tinuum parameters inferred from the RGS data are also 
consistent with those derived from EPIC-PN and EPIC- 
MOS. The intrinsic isotropic luminosities in the 0.3 - 2 
keV and 2-10 keV regions axe Lx ~ 2 x 10 44 erg s 1 
and ~ 5 x 10 43 erg s _1 (assuming H 0 = 65 km s -1 Mpc -1 
and go = 0.5), respectively, which are lower than both the 
ROSAT PSPC and ASCA values by a factor of — 5. The 
measured ion column densities are listed in Table 1. To 
illustrate the statistical significance of the various compo- 
nents, we also list the changes in x 2 when each of the com- 
ponents is removed from the model. The best-fit FWHM 
turbulent velocity is u tU rb = 1430^280 km s J , which is 
much larger than the thermal velocity of a photoion- 
ized medium. We also find weak evidence of an average 
bulk outflow velocity shift with v s hift = +200i}gQ km s 1 , 
where positive velocity denotes a blueshift. 

4. Implications of the Results of Spectral Fitting 

For the measured ion column densities listed in Table 1, 
many of the strong absorption lines of neon and iron L ions 
are in the logarithmic region of the curves of growth, which 
indicates that the derived column densities are highly cou- 
pled with the assumed turbulent velocity. Those of K-shell 
carbon, nitrogen, and oxygen, and iron M-shell ions, how- 
ever, are in the linear regime, and the derived column 
densities are not very sensitive to the turbulent velocity. 

As shown in Table 1, the detections of absorption lines 
from Fe VII - XII and Fexvn - xx are highly significant. 
The column densities of the intermediate charge states 
of Fexm - XVI, however, are consistent with zero. This 
indicates that the line-of-sight material consists of either 
a multi-phase gas in a single medium, or two (or more) 
spatially distinct regions. Motivated by this observational 
fact, we refit the spectrum using the same model, ex- 
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Table 1 . Measured Ion Column Densities 


Ion 

Ni (cm i 

I 

a 

Ax 2 6 

C V 

(6-31^ 

X 

10 1 ® 

4.8 

C VI 

(6.ol®;]) 

X 

10 16 

23.7 

N VI 

(2.41,'q) 

X 

10 16 

29.1 

N VII 

(1.312.5) 

(3-71U) 

X 

10 17 

109.5 

O VII 

X 

10 1 ® 

24.2 

O VIII 

(9.515;;) 

X 

10 1 ® 

24.7 

Ne IX 

(l-2lo,5) 

X 

10 1 ® 

9.9 

NeX 

(4.9i;°,) 

X 

10 17 

25.9 

Fe XVII 

(1-711) 

X 

10 17 

23.0 

Fe XVIII 

(6.3l5;i) 

(9.7l]j) 

X 

10 17 

115.0 

Fe XIX 

X 

10 17 

166.7 

Fe XX 

(3.0+*) 

X 

10 17 

39.5 

Fe VII 

(1 5±i;|) 

X 

10 16 

9.6 

Fe VIII 

(4.6111) 

X 

10 1 ® 

91.8 

FeIX 

(8-8li 2 j) 

(2.411.5) 

X 

10 15 

29.3 

Fe X 

X 

10 1 ® 

72.4 

FeXI 

(i-*iJ:S) 

X 

10 1 ® 

65.3 

Fe XII 

(6.4ll° 4 ) 

X 

10 16 

25.6 


° The turbulent velocities for all of the ions are fixed to the 
same value. We find a best-fit with t>turb = 1 430^280 km s 1 
FWHM. 

6 The increase in \ 2 when the ion component is excluded from 
the best-fit model. 


cept, we assume that the line-of-sight material consists of 
two discrete velocity components; (1) the low-ionization- 
parameter component including Cvi, Cv, Nvi, Ovn, 
and the M-shell Fe, and (2) the high-ionization-parameter 
component including N VII, O VIII, Ne IX, Ne x, and L-shell 
Fe. Contrary to our previous fit where the bulk veloci- 
ties of all the ions were fixed relative to one another, we 
find that the low-ionization-parameter component is sig- 
nificantly blue-shifted with u 8 hift,iow = +420tJgo km s \ 
while the bulk velocity shift in the high-ionization compo- 
nent is consistent with zero (v S hift,high = —20^^ km s~ l ). 
This implies that the low-ionization gas is being acceler- 
ated substantially compared to the high-ionization compo- 
nent, as would be expected in a radiatively-driven outflow 
(Arav & Li 1994). The derived turbulent velocities of both 
components remain consistent with that of our previous 
fit (t>turb,FWHM ~ 1500 km S *). 

From the observed distribution of charge states of M- 
shell iron, the average ionization parameter, £ = L/nr 2 , of 
the absorbing gas is estimated to be log £ ~ 0 based on a 
calculation with the photoionization code XSTAR (Kali- 
man & Krolik 1995) using the inferred continuum shape 
for the illuminating spectrum. The measured ion column 
densities suggest that the corresponding equivalent hydro- 
gen column density is Nh ~ 2 x 10 21 cm" 2 assuming a 
solar iron abundance. This low-£ gas accounts for most 
of the carbon and He-like nitrogen and oxygen absorption 


lines as well. The lack of absorption from Fexm - XVI, 
however, indicates that substantial amounts of material 
with ionization parameters in the range 1 < log £ < 2 are 
not present along the line-of-sight. Whether this is related 
to the global structure of the circumnuclear medium or a 
mere coincidence from a superposition of physically dis- 
tinct regions is not known. 

The absorption lines from H-like nitrogen and oxygen, 
H- and He-like neon, and L-shell iron are produced in a 
medium with 2.0 £ log £ £2.5 and an equivalent hydrogen 
column density of iV H ~ 1 x 10 22 cm -2 . For a normal dust- 
to-gas ratio, the observed reddening of E(B — V) = 0.3 
(Wills et al. 1992) corresponds to a hydrogen column den- 
sity of Nh = 1.7 x 10 21 cm -2 (Burstein & Heiles 1978). 
This value is slightly lower than the total amount of X- 
ray absorbing material observed in the present X-ray spec- 
trum. Coincidentally, however, the derived column density 
of the low-£ component is very close to that of the optical 
reddening, although, we cannot conclusively associate the 
low-£ X-ray absorber with the dusty torus. The column 
density of the high-£ component, on the other hand, is a 
factor of ~ 5 higher. 

An interesting point to note is that the 1 £ log £ ^ 2 
region is not thermally unstable, based on XSTAR cal- 
culations described above. On the other hand, the high-£ 
region (2.0 £ logf <i 2.5) that we observe in the spec- 
trum is thermally unstable. However, complications such 
as non-solar metal abundances and/or inaccurate ioniza- 
tion and recombination rates may alter the shape of the 
thermal stability curve significantly, and, hence, the tem- 
perature ranges of the unstable regions (Hess, Kahn, & 
Paerels 1997; Savin et al. 1999). 


5. Comparisons with Other AGNs 

The absorption spectrum of IRAS 13349+2438 is qualita- 
tively similar to those obtained with the Chandra trans- 
mission grating observations of the Seyfert 1 galaxies 
NGC 5548 and NGC 3783, which show narrow reso- 
nance absorption lines blue-shifted by several hundred 
km s™* 1 . (Kaastra et al. 2000; Kaspi et al. 2000). The de- 
rived ion column densities in these sources, as well as in 
IRAS 13349+2438, are in the range TV* ~ 10 16 -T0 17 cm" 2 , 
and are not high enough to produce observable absorption 
edges. 

Conceptually, the low-ionization component observed 
in IRAS 13349+2438 is similar to the “lukewarm ab- 
sorber” that explains both the observed optical and X- 
ray attenuation in NGC 3227 (Kraemer et al. 2000). 
The spectroscopic signatures, however, are very different. 
In particular, for the column densities derived from the 
IRAS 13349+2438 data, the absorption features are dom- 
inated by discrete resonance line transitions, mainly in 
He-like ions and M-shell iron, and not by photoelectric 
edges as in the model of Kraemer et al. (2000). 
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* • 


As demonstrated in our detailed spectral analysis of 
IRAS 13349+2438, the UTA feature is potentially a pow- 
erful diagnostic tool for probing cool absorbing mate- 
rial using high-resolution X-ray observations. If the low- 
ionization component is indeed associated with the dusty 
torus as the derived column density suggests, this feature 
should be detectable in other AGNs where the line-of-sight 
is partially obscured by the torus. 
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